A within- and across-country assessment of the genomic diversity and autozygosity of South African and eSwatini Nguni cattle.

In southern Africa, the Nguni cattle breed is classified as an indigenous and transboundary animal genetic resource (AnGR) that manifests unique adaptation abilities across distinct agroecological zones. The genetic integrity of various ecotypes is under potential threat due to both indiscriminate crossbreeding and uncontrolled inbreeding. The aim of this study was to assess the genetic diversity and autozygosity that exist both across countries (ES: eSwatini; SA: South Africa) and within countries (SA), between purebred stud animals (SA-S) and research herds (SA-R). Subsets of 96 ES, 96 SA-S, and 96 SA-R genotyped for 40,930 common SNPs were used to study genome-wide profiles of runs of homozygosity (ROH) and heterozygosity (ROHet) as well as inbreeding levels and population structure. The highest percentage (39.8%) of the 2168 ROH segments was 4-8 Mbp in length, whereas 65% of the 935 ROHet segments fell within the 0.5-1 Mbp length category. Inbreeding coefficients indicated positive but low inbreeding (FROH>1Mbp range: 0.025 for SA-S to 0.029 for SA-R). Principal component (PCA) and population structure analyses illustrated genome-level distinctness of (1) the Nguni from global indicine (Boran) and taurine (Hereford) breeds (K = 3), (2) the SA Nguni populations from the ES Nguni population (K = 4), and (3) different Nguni ecotypes within countries (K = 8). Furthermore, greater admixture was observed for the SA-R population compared to purebred SA-S population (shared ancestry = 0.631 ± 0.353 compared to 0.741 ± 0.123), and fewer genomics-defined ES ecotypes were observed than phenotypically (pre)defined. Overall, the results illustrated that genetic uniqueness within the sampled Nguni cattle resulted from both geographic isolation and exposure to different breeding strategies (and, selection pressures). A further loss of genetic variability should be monitored to prevent the endangerment of unique and beneficial ecotypes.

Conservation status and historical relatedness of South African communal indigenous goat populations using a genome-wide single-nucleotide polymorphism marker.

Indigenous goats form the majority of populations in smallholder, low input, low output production systems and are considered an important genetic resource due to their adaptability to different production environments and support of communal farming. Effective population size (N e ), inbreeding levels, and the runs of homozygosity (ROHs) are effective tools for exploring the genetic diversity and understanding the demographic history in efforts to support breeding strategies to use and conserve genetic resources. Across populations, the current N e of Gauteng was the lowest at 371 animals, while the historical N e across populations suggests that the ancestor N e has decreased by 53.86%, 44.58%, 42.16%, and 41.16% in Free State (FS), North West (NW), Limpopo (LP), and Gauteng (GP), respectively, over the last 971 generations. Genomic inbreeding levels related to ancient kinship (F ROH > 5 Mb) were highest in FS (0.08 ± 0.09) and lowest in the Eastern Cape (EC) (0.02 ± 0.02). A total of 871 ROH island regions which include important environmental adaptation and hermo-tolerance genes such as IL10RB, IL23A, FGF9, IGF1, EGR1, MTOR, and MAPK3 were identified (occurring in over 20% of the samples) in FS (n = 37), GP (n = 42), and NW (n = 2) populations only. The mean length of ROH across populations was 7.76 Mb and ranged from 1.61 Mb in KwaZulu-Natal (KZN) to 98.05 Mb (GP and NW). The distribution of ROH according to their size showed that the majority (n = 1949) of the detected ROH were > 5 Mb in length compared to the other categories. Assuming two hypothetical ancestral populations, the populations from KZN and LP are revealed, supporting PC 1. The genomes of KZN and LP share a common origin but have substantial admixture from the EC and NW populations. The findings revealed that the occurrence of high N e and autozygosity varied largely across breeds in communal indigenous goat populations at recent and ancient events when a genome-wide single-nucleotide polymorphism (SNP) marker was used. The use of Illumina goat SNP50K BeadChip shows that there was a migration route of communal indigenous goat populations from the northern part (LP) of South Africa to the eastern areas of the KZN that confirmed their historical relatedness and coincides with the migration periods of the Bantu nation.

Selection signature analysis and genome-wide divergence of South African Merino breeds from their founders.

Merino sheep are a breed of choice across the world, popularly kept for their wool and mutton value. They are often reared as a pure breed or used in crossbreeding and are a common component in synthetic breed development. This study evaluated genetic diversity, population structure, and breed divergence in 279 animals of Merino and Merino-based sheep breeds in South Africa using the Illumina Ovine SNP 50K BeadChip. The sheep breeds analysed included the three Merino-derived breeds of Dohne Merino (n = 50); Meatmaster (n = 47); and Afrino (n = 52) and five presumed ancestral populations of Merinos (Merino (n = 46); South African Merino (n = 10); and South African Mutton Merino (n = 8)); and the non-Merino founding breeds of Damara (n = 20); Ronderib Afrikaner (n = 17); and Nguni (n = 29). Highest genetic diversity values were observed in the Dohne Merino (DM), with H o = 0.39 ± 0.01, followed by the Meatmaster and South African Merino (SAM), with H o = 0.37 ± 0.03. The level of inbreeding ranged from 0.0 ± 0.02 (DM) to 0.27 ± 0.05 (Nguni). Analysis of molecular variance (AMOVA) showed high within-population variance (>80%) across all population categories. The first principal component (PC1) separated the Merino, South African Mutton Merino (SAMM), DM, and Afrino (AFR) from the Meatmaster, Damara, Nguni, and Ronderib Afrikaner (RDA). PC2 aligned each Merino-derived breed with its presumed ancestors and separated the SAMM from the Merino and SAM. The iHS analysis yielded selection sweeps across the AFR (12 sweeps), Meatmaster (four sweeps), and DM (29 sweeps). Hair/wool trait genes such as FGF12; metabolic genes of ICA1, NXPH1, and GPR171; and immune response genes of IL22, IL26, IFNAR1, and IL10RB were reported. Other genes include HMGA, which was observed as selection signatures in other populations; WNT5A, important in the development of the skeleton and mammary glands; ANTXR2, associated with adaptation to variation in climatic conditions; and BMP2, which has been reported as strongly selected in both fat-tailed and thin-tailed sheep. The DM vs. SAMM shared all six sweep regions on chromosomes 1, 10, and 11 with AFR vs. SAMM. Genes such as FGF12 on OAR 1:191.3-194.7 Mb and MAP2K4 on OAR 11:28.6-31.3 Mb were observed. The selection sweep on chromosome 10 region 28.6-30.3 Mb harbouring the RXFP2 for polledness was shared between the DM vs. Merino, the Meatmaster vs. Merino, and the Meatmaster vs. Nguni. The DM vs. Merino and the Meatmaster vs. Merino also shared an Rsb-based selection sweep on chromosome 1 region 268.5-269.9 Mb associated with the Calpain gene, CAPN7. The study demonstrated some genetic similarities between the Merino and Merino-derived breeds emanating from common founding populations and some divergence driven by breed-specific selection goals. Overall, information regarding the evolution of these composite breeds from their founding population will guide future breed improvement programs and management and conservation efforts.

Bacterial profiling of Haemonchus contortus gut microbiome infecting Dohne Merino sheep in South Africa.

A metagenomic approach was used to study the gut microbiome of Haemonchus contortus field strains and that of its predilection site, the abomasum of Dohne Merino sheep. The abomasum contents and H. contortus were collected from 10 naturally infected Dohne Merino sheep. The H. contortus specimens were classified and sexually differentiated using morphometric characters and was further confirmed through molecular identification. We investigated differences and similarities between the bacterial composition of the adult male and female H. contortus gut microbiomes, which were both dominated by bacteria from the Escherichia, Shigella, Vibrio and Halomonas genera. Major abundance variations were identified between the shared adult male and female H. contortus microbiomes. The results also revealed that Succiniclasticum, Rikenellaceae RC9 gut group and Candidatus Saccharimonas were the predominant genera in the Dohne Merino abomasum. This study provides insight into the highly diverse bacterial composition of the H. contortus gut microbiome and the Dohne Merino abomasum which needs to be studied further to explore the complex interactions of different gastrointestinal nematode microbiomes with the host.

Runs of homozygosity analysis of South African sheep breeds from various production systems investigated using OvineSNP50k data.

BACKGROUND: Population history, production system and within-breed selection pressure impacts the genome architecture resulting in reduced genetic diversity and increased frequency of runs of homozygosity islands. This study tested the hypothesis that production systems geared towards specific traits of importance or natural or artificial selection pressures influenced the occurrence and distribution of runs of homozygosity (ROH) in the South African sheep population. The Illumina OvineSNP50 BeadChip was used to genotype 400 sheep belonging to 13 breeds from South Africa representing mutton, pelt and mutton and wool dual-purpose breeds, including indigenous non-descript breeds that are reared by smallholder farmers. To get more insight into the autozygosity and distribution of ROH islands of South African breeds relative to global populations, 623 genotypes of sheep from worldwide populations were included in the analysis. Runs of homozygosity were computed at cut-offs of 1-6 Mb, 6-12 Mb, 12-24 Mb, 24-48 Mb and > 48 Mb, using the R package detectRUNS. The Golden Helix SVS program was used to investigate the ROH islands. RESULTS: A total of 121,399 ROH with mean number of ROH per animal per breed ranging from 800 (African White Dorper) to 15,097 (Australian Poll Dorset) were obtained. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were in the short (1-6 Mb) category (88.2%). Most animals had no ROH > 48 Mb. Of the South African breeds, the Nguni and the Blackhead Persian displayed high ROH based inbreeding (FROH) of 0.31 ± 0.05 and 0.31 ± 0.04, respectively. Highest incidence of common runs per SNP across breeds was observed on chromosome 10 with over 250 incidences of common ROHs. Mean proportion of SNPs per breed per ROH island ranged from 0.02 ± 0.15 (island ROH224 on chromosome 23) to 0.13 ± 0.29 (island ROH175 on chromosome 15). Seventeen (17) of the islands had SNPs observed in single populations (unique ROH islands). The MacArthur Merino (MCM) population had five unique ROH islands followed by Blackhead Persian and Nguni with three each whilst the South African Mutton Merino, SA Merino, White Vital Swakara, Karakul, Dorset Horn and Chinese Merino each had one unique ROH island. Genes within ROH islands were associated with predominantly metabolic and immune response traits and predomestic selection for traits such as presence or absence of horns. CONCLUSIONS: Overall, the frequency and patterns of distribution of ROH observed in this study corresponds to the breed history and implied selection pressures exposed to the sheep populations under study.

The genomic architecture of South African mutton, pelt, dual-purpose and nondescript sheep breeds relative to global sheep populations.

South Africa has a diverse array of phenotypically distinct and locally adapted sheep breeds that have been developed for different production systems ranging from mutton to wool and pelt, and some dual-purpose and nondescript breeds kept by smallholder farmers. This study investigated genetic diversity, population genetic structure and divergence between South African sheep breeds in order to gain an insight into breed history and genomic architecture aligned to breeding goals and production systems. The Illumina OvineSNP50 BeadChip was used to genotype 400 sheep belonging to 14 breeds representing mutton, pelt and mutton and wool dual-purpose breeds. Nguni sheep were included as a representative of indigenous nondescript breeds that are reared by smallholder farmers. Seeking a clearer understanding of the genetic diversity of South African breeds relative to global populations, 623 genotypes of sheep from worldwide populations were included in the analysis. These sheep breeds included six African, two Asian and eight European breeds. Across breeds, genetic diversity ranged from observed heterozygosity (H0 ) = 0.26 ± 0.02 in Namaqua Afrikaner to H0  = 0.38 ± 0.01 in Dohne Merino. The overall mean H0 was 0.35 ± 0.04. The African and Asian populations were the most inbred populations with FIS ranging from 0.17 ± 0.05 in Grey Swakara and Ronderib Afrikaner sheep to 0.34 ± 0.07 in the Namaqua Afrikaner. The South African Dohne Merino (FIS  = 0.03 ± 0.01), SA Merino (FIS  = 0.05 ± 0.04) and Afrino (FIS  = 0.09 ± 0.02) and other global Merino-derived breeds were the least inbred. The first principal component explained 27.7% of the variation and separated the fat- and rump-tailed sheep (i.e. Swakara, Nguni, Blackhead Persian, Ethiopian Menzi, Meatmaster) from the Merino and Merino-derived breeds and the Dorset Horn. The second principal component separated the Merino and Merino-derived breeds from the English breed of Dorset Horn. Overall, South African indigenous breeds clustered together with indigenous breeds from other African and Asian countries. The optimal admixture cluster (K = 20) revealed various sources of within- and amongst-breed genomic variation associated with production purpose, adaptation and history of the breeds. The Blackhead Persian, Nguni and Namaqua Afrikaner breeds differed significantly from other breeds, particularly with the South African Mutton Merino and Dorset Horn. Breed-differentiating SNPs were observed within genomic regions associated with growth, adaptation and reproduction. Genes such as RAB44, associated with growth and meat/carcass traits, differentiated the Blackhead Persian from the Dorset Horn and South African Mutton Merino. The MAP2 and HRAS genes, which are associated with immune traits involving Toll-like receptors and Chemokine signalling pathways, differentiated the Nguni from the Dorset Horn. The current results give insight into the current status of the sheep genetic resources of South Africa relative to the global sheep population, highlighting both genetic similarities as well as divergence associated with production systems, geographical distribution and local adaptation.

Population structure of indigenous southern African goats based on the Illumina Goat50K SNP panel.

In this study, the genetic structure of indigenous Tswana and Swazi goats using the Illumina Goat50K SNP array was investigated. Two South African commercial goat breeds were included to investigate admixture with the indigenous populations in southern Africa. A total of 144 DNA samples including Boer goats (n = 24), Kalahari Red (n = 24), Swazi (n = 48), and Tswana goats (n = 48) were genotyped. Statistical analysis was performed using PLINK version 1.07. Genetic diversity, measured as expected heterozygosity, was estimated at 0.390, 0.398, 0.413, and 0.387 for Boer, Kalahari Red, Tswana, and Swazi goats, respectively. The individual inbreeding coefficient varied from 0.019 ± 0.05 to 0.011 ± 0.06 for the Tswana and Swazi goats, respectively. The Principal component analysis clustered the populations according to geographical origin and breed type. Linkage disequilibrium (LD) for shorter intervals (0-10 kb) ranged from 0.44 to 0.56 and commercial breeds had higher values. Effective population sizes decreased with generations and at the 13th generation ranged between 87 for Boer to 266 for Tswana goats. The Tswana population exhibited the highest level of genetic variation and effective population size, which holds potential for improved production in marginal regions. A national strategy is required to maintain genetic diversity in communal goat production systems through well-structured breeding and conservation programs.

Relationship between population genomic structure and growth profiles of South African goats under different production systems.

Goats play a major role in poor marginalized communities of South Africa for food security and socio-economic purposes. Majority of the goats are raised in villages with poor infrastructure and resources, therefore facing challenges that affect growth performance which leads to low mature weights. Investigating growth profiles will shed light on growth performances and will aid in goat improvement and selection. This study investigated the growth profiles and genomic structure of SA indigenous breeds raised in different production systems to unravel the genetic potential of indigenous goat populations. Live weights and morphological body measurements were collected from a total of 83 kids representing the commercial meat-producing SA Boer (n = 14); the indigenous veld goats (IVG) of NC Skilder (n = 14), Mbuzi (n = 13), and Xhosa lob (n = 14) raised under intensive systems; and nondescript village goat populations (n = 14) raised in intensive and others (n = 14) raised in extensive production systems. The remaining 72 of 83 phenotyped goats were genotyped using the Illumina Caprine SNP50K BeadChip. The SA Boer had a higher weight (28.96 ± 0.30 kg) gain as compared to other populations. The Mbuzi population was the smallest (14.83 ± 0.33 kg), while the village goats raised in Pella Village were relatively smaller (17.55 ± 0.37 kg) than those raised on the research farm (19.55 ± 0.36 kg). The study concluded that both genetics and management systems can lead to improved growth performance in goat production. The outputs of this study can be used to identify suitable breeds and potential genotypes for optimal growth and establish optimal goat management systems suitable for communal farmers for improved productivity.

Gene expression profiles of the small intestine of village chickens from an Ascaridia galli infested environment.

Nematodes of the genus Ascaridia are known to infect many species of birds and result in fatal diseases. A. galli damages the intestinal mucosa of chickens leading to blood loss, secondary infection and occasionally the obstruction of small intestines due to high worm burden. This study investigated the gene expression profiles in chickens from two different provinces of South Africa naturally exposed to A. galli infestations and tested either positive or negative for the parasite. The study further investigated gene expression profiles of the A. galli infected duodenum, jejunum and ileum tissues of the small intestines. The A. galli positive intestines displayed hypertrophy of the intestinal villi with accumulation of inflammatory cells and necrosis of the crypts of Lieberkühn glands, lesions that were absent in the uninfected intestines. Total RNA isolated from small intestines of infected and non-infected intestines was sequenced using Illumina HiSeq technology to generate up to 23,856,130 reads. Between any two-way comparisons of the intestines, 277 and 190 transcripts were significantly expressed in Limpopo and KwaZulu-Natal (KZN) chickens, respectively. Gene ontology analysis of the differentially expressed genes (DEGs) revealed an enrichment of genes reported to function in the immune response, defense response, inflammatory response and cell signalling genes. T cell receptor signalling pathways and arachidonic acid metabolism pathways were among the most significantly impacted pathways. Overall, the study provided insights into adaptative mechanisms for chickens extensively raised in parasite infected environments.

Gene expression profiles of the small intestine of village chickens from an Ascaridia galli infested environment.

Nematodes of the genus Ascaridia are known to infect many species of birds and result in fatal diseases. A. galli damages the intestinal mucosa of chickens leading to blood loss, secondary infection and occasionally the obstruction of small intestines due to high worm burden. This study investigated the gene expression profiles in chickens from two different provinces of South Africa naturally exposed to A. galli infestations and tested either positive or negative for the parasite. The study further investigated gene expression profiles of the A. galli infected duodenum, jejunum and ileum tissues of the small intestines. The A. galli positive intestines displayed hypertrophy of the intestinal villi with accumulation of inflammatory cells and necrosis of the crypts of Lieberkühn glands, lesions that were absent in the uninfected intestines. Total RNA isolated from small intestines of infected and non-infected intestines was sequenced using Illumina HiSeq technology to generate up to 23,856,130 reads. Between any two-way comparisons of the intestines, 277 and 190 transcripts were significantly expressed in Limpopo and KwaZulu-Natal (KZN) chickens, respectively. Gene ontology analysis of the differentially expressed genes (DEGs) revealed an enrichment of genes reported to function in the immune response, defense response, inflammatory response and cell signalling genes. T cell receptor signalling pathways and arachidonic acid metabolism pathways were among the most significantly impacted pathways. Overall, the study provided insights into adaptative mechanisms for chickens extensively raised in parasite infected environments.

Landscape genomics and pathway analysis to understand genetic adaptation of South African indigenous goat populations.

In Africa, extensively raised livestock populations in most smallholder farming communities are exposed to harsh and heterogeneous climatic conditions and disease pathogens that they adapt to in order to survive. Majority of these livestock species, including goats, are of non-descript and uncharacterized breeds and their response to natural selection presented by heterogeneous environments is still unresolved. This study investigated genetic diversity and its association with environmental and geographic conditions in 194 South African indigenous goats from different geographic locations genotyped on the Illumina goat SNP50K panel. Population structure analysis revealed a homogeneous genetic cluster of the Tankwa goats, restricted to the Northern Cape province. Overall, the Boer, Kalahari Red, and Savanna showed a wide geographic spread of shared genetic components, whereas the village ecotypes revealed a longitudinal distribution. The relative importance of environmental factors on genetic variation of goat populations was assessed using redundancy analysis (RDA). Climatic and geographic variables explained 22% of the total variation while climatic variables alone accounted for 17% of the diversity. Geographic variables solitarily explained 1% of the total variation. The first axis (Model I) of the RDA analysis revealed 329 outlier SNPs. Landscape genomic approaches of spatial analysis method (SAM) identified a total of 843 (1.75%) SNPs, while latent factor mixed models (LFMM) identified 714 (1.48%) SNPs significantly associated with environmental variables. Significant markers were within genes involved in biological functions potentially important for environmental adaptation. Overall, the study suggested environmental factors to have some effect in shaping the genetic variation of South African indigenous goat populations. Loci observed to be significant and under selection may be responsible for the adaption of the goat populations to local production systems.

Management and control of gastrointestinal nematodes in communal goat farms in Zimbabwe.

Goats are an important source of livelihood especially in smallholder communities. Infections with gastrointestinal nematodes (GIN) remain the most prevalent parasitic diseases affecting small ruminants. The study was conducted to assess management, the level of knowledge and control of gastrointestinal nematodes. Surveys were conducted in Chipinge, Shurugwi, Binga, Tsholotsho and Matobo districts, representing the five natural/agro-ecological regions (NR) in Zimbabwe. Data was collected in 135 households using a pre-tested semi-structured questionnaire. Results indicated that goats were ranked the most important livestock species, with high flock sizes in NR IV and V. Partitioning of roles was such that the adult males were involved in decision-making while females and children were involved in day-to-day management of animals. Farmers showed low levels of input use, with natural pasture (98.4%) being the main feed source and indigenous breeds (73.2%) being kept. Farmers ranked food and financial benefits as the main reasons for keeping goats. Gastrointestinal nematodes ranked the highest as the most common disease, with majority of farmers (57%) not controlling or treating animals and 63% of farmers not having knowledge on the spread of GIN. Access to veterinary services, anthelmintic class used and breeds used by the farmers had the highest effects on parasitic infections in households. Farmer education is required for capacitation of farmer in terms of disease prevention and control so as to improve goat production.

Prevalence and risk factors of gastrointestinal parasitic infections in goats in low-input low-output farming systems in Zimbabwe.

A longitudinal study was conducted in low-input low-output farming systems to determine the prevalence of gastrointestinal parasitic infections in different age groups, sex and associated risk factors in goats. A total of 580 indigenous goats were randomly selected in areas representing the five agro-ecological regions of Zimbabwe in the dry and wet seasons. Blood and faecal samples were collected from each animal and egg/oocyst per gram of faeces (epg/opg), larval culture, and packed cell volumes (PCV) were determined. Factors affecting parasitic infections were evaluated. Highest prevalence was determined for Eimeria oocysts (43%), strongyles (31%) and lower levels in trematodes and cestodes. Parasites identified were Haemonchus, Strongyloides and Oesophagostomum. Area, season, sex and age significantly influenced patterns of gastrointestinal infections (P < 0.05). Cannonical correlations indicated that parasite species composition varied by area and impacts of risk factors also differed. Risk of infection was very high for goats sampled in Natural regions (NR) I, II, III (OR = 6.6-8.2; P < 0.05) as compared to those in NR IV and V. Highest helminths and Eimeria infections were observed in the wet vs. dry season (P < 0.05). Young animals were more susceptible to parasitic infections (P < 0.05). Prevalence was higher in males than females, with odds of infection for males being almost three times to that for females (P < 0.0001). Knowledge concerning gastrointestinal helminth biology and epidemiological infection patterns caused by these parasites is essential in the development of appropriate control strategies and this has a potential to reduce production losses.

Breeding for resistance to gastrointestinal nematodes - the potential in low-input/output small ruminant production systems.

The control of gastrointestinal nematodes (GIN) is mainly based on the use of drugs, grazing management, use of copper oxide wire particles and bioactive forages. Resistance to anthelmintic drugs in small ruminants is documented worldwide. Host genetic resistance to parasites, has been increasingly used as a complementary control strategy, along with the conventional intervention methods mentioned above. Genetic diversity in resistance to GIN has been well studied in experimental and commercial flocks in temperate climates and more developed economies. However, there are very few report outputs from the more extensive low-input/output smallholder systems in developing and emerging countries. Furthermore, results on quantitative trait loci (QTL) associated with nematode resistance from various studies have not always been consistent, mainly due to the different nematodes studied, different host breeds, ages, climates, natural infections versus artificial challenges, infection level at sampling periods, among others. The increasing use of genetic markers (Single Nucleotide Polymorphisms, SNPs) in GWAS or the use of whole genome sequence data and a plethora of analytic methods offer the potential to identify loci or regions associated nematode resistance. Genomic selection as a genome-wide level method overcomes the need to identify candidate genes. Benefits in genomic selection are now being realised in dairy cattle and sheep under commercial settings in the more advanced countries. However, despite the commercial benefits of using these tools, there are practical problems associated with incorporating the use of marker-assisted selection or genomic selection in low-input/output smallholder farming systems breeding schemes. Unlike anthelmintic resistance, there is no empirical evidence suggesting that nematodes will evolve rapidly in response to resistant hosts. The strategy of nematode control has evolved to a more practical manipulation of host-parasite equilibrium in grazing systems by implementation of various strategies, in which improvement of genetic resistance of small ruminant should be included. Therefore, selection for resistant hosts can be considered as one of the sustainable control strategy, although it will be most effective when used to complement other control strategies such as grazing management and improving efficiency of anthelmintics currently.

Population genomic structure and linkage disequilibrium analysis of South African goat breeds using genome-wide SNP data.

The sustainability of goat farming in marginal areas of southern Africa depends on local breeds that are adapted to specific agro-ecological conditions. Unimproved non-descript goats are the main genetic resources used for the development of commercial meat-type breeds of South Africa. Little is known about genetic diversity and the genetics of adaptation of these indigenous goat populations. This study investigated the genetic diversity, population structure and breed relations, linkage disequilibrium, effective population size and persistence of gametic phase in goat populations of South Africa. Three locally developed meat-type breeds of the Boer (n = 33), Savanna (n = 31), Kalahari Red (n = 40), a feral breed of Tankwa (n = 25) and unimproved non-descript village ecotypes (n = 110) from four goat-producing provinces of the Eastern Cape, KwaZulu-Natal, Limpopo and North West were assessed using the Illumina Goat 50K SNP Bead Chip assay. The proportion of SNPs with minor allele frequencies >0.05 ranged from 84.22% in the Tankwa to 97.58% in the Xhosa ecotype, with a mean of 0.32 ± 0.13 across populations. Principal components analysis, admixture and pairwise FST identified Tankwa as a genetically distinct population and supported clustering of the populations according to their historical origins. Genome-wide FST identified 101 markers potentially under positive selection in the Tankwa. Average linkage disequilibrium was highest in the Tankwa (r(2)  = 0.25 ± 0.26) and lowest in the village ecotypes (r(2) range = 0.09 ± 0.12 to 0.11 ± 0.14). We observed an effective population size of <150 for all populations 13 generations ago. The estimated correlations for all breed pairs were lower than 0.80 at marker distances >100 kb with the exception of those in Savanna and Tswana populations. This study highlights the high level of genetic diversity in South African indigenous goats as well as the utility of the genome-wide SNP marker panels in genetic studies of these populations.

Targeted high-throughput growth hormone 1 gene sequencing reveals high within-breed genetic diversity in South African goats.

This study assessed the genetic diversity in the growth hormone 1 gene (GH1) within and between South African goat breeds. Polymerase chain reaction-targeted gene amplification together with Illumina MiSeq next-generation sequencing (NGS) was used to generate the full length (2.54 kb) of the growth hormone 1 gene and screen for SNPs in the South African Boer (SAB) (n = 17), Tankwa (n = 15) and South African village (n = 35) goat populations. A range of 27-58 SNPs per population were observed. Mutations resulting in amino acid changes were observed at exons 2 and 5. Higher within-breed diversity of 97.37% was observed within the population category consisting of SA village ecotypes and the Tankwa goats. Highest pairwise FST values ranging from 0.148 to 0.356 were observed between the SAB and both the South African village and Tankwa feral goat populations. Phylogenetic analysis indicated nine genetic clusters, which reflected close relationships between the South African populations and the other international breeds with the exception of the Italian Sarda breeds. Results imply greater potential for within-population selection programs, particularly with SA village goats.

Population genetic structure of Ascaridia galli of extensively raised chickens of South Africa.

Ascaridia galli is one of the most common nematode affecting chickens. This study characterized A. galli parasites collected from South African village chickens of Limpopo (n=18) and KwaZulu-Natal (n=22) provinces using the 510bp sequences of cytochrome C oxidase subunit 1 gene of the mitochondrial DNA. Fourteen and 12 polymorphic sites were observed for Limpopo and KwaZulu-Natal sequences, respectively. Six haplotypes were observed in total. Haplotype diversity was high and ranged from 0.749 for Limpopo province to 0.758 for KwaZulu-Natal province isolates. There was no genetic differentiation between A. galli from Limpopo and KwaZulu-Natal provinces. The six South African haplotypes were unique compared to those published in the GeneBank sampled from Hy-line chickens raised under organic farming in Denmark. The utility of cytochrome C oxidase subunit 1 gene as a potential genetic marker for studying A. galli in village chicken populations is presented.

Tropically adapted cattle of Africa: perspectives on potential role of copy number variations.

Africa is host to diverse and locally adapted cattle breeds that are expected to survive the harsh and extreme tropical environments associated with diseases and parasite infections, heat stress and episodes of feed and water scarcity. Genomic copy number variations (CNVs) are considered to be primary role players in cattle breed formation and adaptation where isolation and genetic drift together with subsequent mutations have created an enormous diversity of local populations. CNVs are modifications in DNA structure comprising deletions, duplications and insertions that are >1 kb in size. Despite attracting much attention, the frequency and pattern of bovine CNV events, especially in African cattle breeds, are for the most part largely unknown. Characterization of genetic variation in the indigenous cattle of Africa will be a vital step toward dissecting the molecular mechanisms underlying phenotypic variation and local adaptation. This review therefore aims to describe the current knowledge regarding bovine CNVs and the implications and potentials they encompass for dissecting genetic adaptation and the genotypic skeleton of tropical African cattle populations.

Haplo-block structure of Southern African village chicken populations inferred using genome-wide SNP data.

This study investigated the haplo-block structure, haplotype sharing, and diversity in extensively raised chicken populations of Southern Africa. Two hundred ninety village chickens from Malawi (N = 30), South Africa (N = 132), and Zimbabwe (N = 128) were included in the study, from which 649, 2104, and 2442 haplo-blocks were observed, respectively. The majority of haplo-blocks were smaller than 25 kb in size and only five blocks were more than 2000 kb in size. The low chromosomal coverage of haplo-blocks observed across the genome suggests that multiple recombination events fragmented the ancestral haplo-blocks into smaller sizes. Haplo-block sharing was observed between populations with 2325 haplo-blocks common between Zimbabwe and Malawi and 2689 between South Africa and Zimbabwe. Haplotype sharing allows transferability of genomic tools between these extensively raised chicken populations of Southern Africa. The unique haplo-blocks could have originated from isolated evolution taking place in specific agro-ecological zones. Quantitative trait loci analysis revealed that genes related to body composition were spanned by these haplo-blocks. Body composition traits are important for village chicken populations, which have to harness poor quality feed obtained from the environment to meet their maintenance and production needs.

Diversity and origin of South African chickens.

The objectives of this study were to analyze the genetic diversity and structure of South African conserved and field chicken populations and to investigate the maternal lineages of these chicken populations. Four South African conserved chicken populations (n = 89), namely, Venda (VD_C), Ovambo, Naked Neck, and Potchefstroom Koekoek from the Animal Production Institute of the Agricultural Research Council, and 2 field populations, the Venda and Ovambo (OV_F), from which the Ovambo and the Venda conservation flocks were assumed to have been sampled, were genotyped for 460 bp of the mitochondrial DNA (mtDNA) D-loop sequence. Haplotypes of these chickens were aligned to 7 Japanese and 9 Chinese and Eurasian chicken mtDNA D-loop sequences taken from GenBank and reflecting populations from presumed centers of domestication. Sequence analysis revealed 48 polymorphic sites that defined 13 haplotypes in the South African chicken populations. All 6 South African conserved and field chicken populations observed were found to be polymorphic, with the number of haplotypes ranging from 3 for VD_C to 8 for OV_F. The lowest haplotype diversity, 0.54 ± 0.08, was observed in VD_C chickens, whereas the highest value, 0.88 ± 0.05, was observed in OV_F chickens. Genetic diversity between the 4 South African conserved and 2 field chicken populations constituted 12.34% of the total genetic variation, whereas within-population diversity constituted 87.66% of the total variation. The median network analysis of the mtDNA D-loop haplotypes observed in the South African conserved and field populations and the reference set resulted in 5 main clades. All 6 South African chickens were equally represented in the major clade, E, which is presumed to be of Indian subcontinent maternal origin and may have its roots in Southeast Asia. The results showed multiple maternal lineages of South African chickens. Conservation flocks and field chicken populations shared the major haplotypes A, D and E, which were presumed to be of Chinese, Southeast Asian, and Indian subcontinental origin.